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Abstract. 

We report on our recent evaluation of the s-channel cut contribution to J /if) hadro- 
production. We show that it is likely significantly larger than the usual cut contribution 
of the colour-singlet model (CSM), which is known to underestimate the experimental 
measurements. Here the s-channcl cut develops for configurations with off-shell quarks in 
the bound state. A correct treatment of its contribution requires the introduction of a four- 
point function, partially constrained by gauge invariance and limiting behaviours at small 
and large momenta. When the unconstrained degrees of freedom are fixed to reproduce the 
Tevatron data, we show that RHIC data are remarkably well reproduced down to very low 
transverse momenta Pt without need of resummation of initial-state gluon effects. This 
unique feature might be typical of s-channel cut contribution. 
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1. INTRODUCTION 

More than ten years after the first measurements by the CDF Collaboration of 
the direct production of J /if) and if)' at y/s = 1.8 TeV [U [2] we are still facing 
disagreements between theoretical predictions from the various available models 
and experimental studies of the cross section and the polarisation from the Tevatron 
and RHIC (for reviews see [3]). 

CDF [3] recently confirmed their previous polarisation measurement [5] showing 
an unpolarised or slightly longitudinally polarised prompt J /if) yield. This has 
reinforced the doubts cast on the dominance of the Colour Octet Mechanism 
(COM) coming from the application of NRQCD |6J. At the same time, many 
new results became available, e.g. the long-awaited NLO QCD corrections to the 
CSM [7] -showing significant enhancement of the cross section, an up-to-date 
proof [8] of NRQCD factorisation; an improved treatment of NRQCD factorisation 
in fragmentation regions where three heavy quarks have similar momenta jU] and 
last but not least a recent evaluation of the dominant a| (NNLO) [10] correction 
to T(nS) production, the latter solving the longstanding conflict between the 
experimental measurements from the Tevatron at mid and large Pt [TlT [12] and 
the prediction from the CSM [13]. 

Considering that none of the existing theoretical approaches could reproduce all 



available experimental data, we undertook in |14J a systematic study of the cut 
contributions due to off-shell and non-static quarks. In particular, we questioned 
the assumption of the CSM that takes the heavy quarks forming the quarkonium 
(<2) as being on-shell p2]. If they are not, the usual s-channel cut contributes to 
the imaginary part of the amplitude and need to be considered on the same footing 
as the CSM cut. 

Current conservation for such off-shell configuration responsible for the s-channel 
cut imposes the introduction of an additional four-point function, or contact 
current [15], accounting for the interactions between the cc pair emitting the 
external gluon. In fact, this mechanism arises because of the possibility that the 
outgoing gluon is emitted by the particle interacting in the dressed c — c — J/ip 
vertex (see Figure [I] (a)), as depicted in Figure [l] (b). The pair of on-shell quarks 
that makes the final J/^-gluon state is now in a colour-octet state which thus 
recovers the necessity for such configurations as a natural consequence of restoring 
gauge invar iance. 

Although current conservation imposes the introduction of such a 4-point vertex, 
it does not enable to relate it univocally to the 3-point one. Yet, there exist certain 
minimal requirements [15] 116] which such a 4-point function should satisfy. The 
4-point function proposed in |T4] provided a conserved current but was not entirely 
satisfactory since it contained poles (by construction similar to the basic direct and 
crossed contributions), and such poles for the contact current are unphysical and 
therefore should be avoided [13] . 

Another caveat to avoid was formerly identified by Drell and Lee [TS]. Indeed, the 
minimal substitution prescription — > d^ + iQA^ (Q: charge, A^: vector potential) 
in an effective Lagrangian corresponding to the dressed hadronic vertex is deficient 
in that it violates the high-energy scaling behaviour, because in avoiding poles for 
the 4-point function, it partially replaces the true momentum dependence of the 
vertices by constants. 

Such an issue can be easily avoided in the approach which we shall follow and 
which was applied to pion photoproduction processes [T71 [T8J [TTJJ [20] . As we shall 
show in the following, it is hence possible to build a 4-point vertex encompassing 
two limiting behaviours, when the final-state gluon is soft or hard [21J. In turn, we 
shall show that this enables to reproduce experimental data from the Tevatron up 
to mid Pt by adjusting the unconstrained parameters of the 4-point vertex and 
hence to get a remarkable agreement with data from RHIC down to low Py. 

2. OUR APPROACH 
2.1. The three-point function 

We shall follow the approach developed in [TJ], where the transition qq — > Q is 
described by the 3-point function 



Tf\p,P) = T(p,P) 1 ,, 



(1) 



where P = p\ — p 2 and p = {p\ + p 2 )/2 are the total and relative momenta, respec- 
tively, of the two quarks bound as a quarkonium state, with p\ and p 2 being their 
individual four-momenta. Ansatz Q amounts to representing the vector meson as 
a massive photon with a non-local coupling. The generic picture of the physical 
origin of the dressed vertex function T(p,P) is given in Figure [jja). In the present 
work, we describe the relative- momentum distribution T(p, P) of the quarks phe- 
nomenologically as a Gaussian form, function of the square of the relative cm. 
3-momentum p of the quarks, which can be written in a Lorentz invariant form as 

p 2 = —p 2 + jrp-. Explicitly, we have 

T(p,P) = Ne-£, (2) 

with a normalisation N -fixed by the leptonic-decay width [E] and a size param- 
eter A, which can be obtained from studies in relativistic quark models [23J. 
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FIGURE 1. Illustration of the mechanisms (a) contributing to the dressing of the 3-point 
function and (b) responsible for the 4-point function I^ 4 ': the external gluon is attached 
here within gluon loops of the dressed vertex, thus producing a diagram without poles and with 
a kinematic behaviour genuinely different from the initial 3-point vertex. 



2.2. The four-point function: minimal substitution 

Instead of directly employing the generalised Ward-Takahashi relations [H] for 
the complete current, we will rather make use of an equivalent condition for the 
contact current similarly to what was done in [T7] in the pion-photoproduction 
case. 

First let us write the 4-point function depicted in Figure [T^b) as 

r« = -ig.T^f , (3) 

where g s is the strong coupling constant, Tf k the colour matrix, and fi and v are 
the Lorentz indices of the outgoing J/ ip and gluon, respectively. For simplicity, we 
have suppressed all indices on the left-hand side. The c — c—J/ip vertex function 
r^ 3 ) with the kinematics of the direct graph is denoted here by T± and for the 
crossed graph by T 2 , i.e., Ti =T (ci — ^ ,P) and r 2 = T (c 2 + f as shown in 
Figures. [3ta) and (b). The gauge-invariance condition for the contact current 



for the outgoing gluon with momentum q reads now 

QuM" = Tx — T 2 (4) 

The contact current can be now constructed as usual [T7J [18], HU ;20j in terms of 
an auxiliary function F = F(c\,C2,q) which contains the remaining unconstrained 
degrees of freedom of the problem. This gives 

(2c2 + gT(IW) (2 Cl - gf (r 2 - F) 
{C2 + q) — m z {Ci — q) z — m z 

where we take c 2 = c 2 = m 2 and P 2 = M 2 from the beginning, with m and M being 
the masses of the quark and the J/ip, respectively. One easily verifies that this 
additional contact current satisfies the gauge-invariance condition Q. 

F(ci,C2,q) has now to be chosen so that the current ^ satisfies crossing sym- 
metry (i.e., symmetry under the exchange c% <-> — c 2 ) and is free of singularities. 
Defining the constant To as the (unphysical) value of the momentum distribution 
T(p,P) when all three legs of the vertex are on their respective mass shells, we 
should then have F(ci,c 2 ,q) = T when (c 2 + g) 2 = m 2 or (ci — q) 2 = m 2 . 

A priori, F = Tq everywhere should be satisfactory. However, this corresponds to 
the minimal substitution discussed by Drell and Lee [IS] and later by Ohta [22] . As 
we mentioned above, this does not provide the correct scaling properties at large 
energies, which means within the present context that F = T would not lead to 
the expected Pt scaling of the amplitude. Numerically, this choice overshoots the 
experimental data by more than one order of magnitude at Pt = 20 GeV as shown 
on Figure [2} 
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FIGURE 2. Comparison between the J/ip production cross section from s-channel cut ob- 
tained with the minimal substitution (F — To) for the 4-point function and the CDF data [2]. 
See below for details on the cross-section evaluation. 



2.3. The four-point function: our proposal 



To avoid such abnormal scaling behaviour, we have to impose, in the large 
relative-momentum region, that the contact term and therefore the function 
F(ci,C2,q) exhibit a fall-off similar to the 3-point vertex functions. 

The simplest crossing-symmetric choice satisfying this property is 



F = r 1 + r 2 -^. (6) 

J- o 

The solution we propose here is to build F(ci,C2,q) from these two limiting cases: 

F =r at low momentum, 

F —Y\ + T 2 — at large momentum. 

To 

To this end, it is practical to choose the following simple Ansatz 

F(c 1 ,c 2 ,g) = r -Mc 1 -c 2 ) (ro - r ^ r °- r2) , (7) 

where the (crossing-symmetric) function /i(ci-c 2 ) rises to become unity for large 
relative momentum. 

Note that F = T at the poles is satisfied independently of h{c\ ■ c 2 ). Indeed, it is 
multiplied on the right by a factor which vanishes at the poles since either Ti = T 
or T 2 = rV A phenomenological choice for the interpolating function h(c\ ■ c 2 ) we 
can then propose is 

K 2 

h( Cl -c 2 ) = l -a— — ■ it, (8) 

K — (ci • C2 + m z ) 

with two parameters, a and k. This choice is in no way unique: in a manner of 
speaking, this choice is simply a way of parameterising our ignorance by employing 
minimal properties of r^. 

Other choices could be analysed but our conclusions, that s-channel cut contri- 
bution can be large and can indeed help to reproduce the experimental data, would 
not be affected. 



3. RESULTS 

For the Tevatron and RHIC kinematics, the direct J/ip are produced by gluon 
fusion and a final-state gluon emission is required to conserve C-parity and provide 
the J/ip with its Pt- The relevant diagrams for the s-channel cut of the LO gluon 
fusion process are shown on Figure |3| We use the same normalisation of as 
in [H], m c = 1.87 GeV and A = 1.8 GeV. As shown in |14j . A can be varied between 
1.2 and 2.2 GeV without affecting much the results. The same statement holds here. 




FIGURE 3. (a) & (b) Leading-order (LO) s-channel cut diagrams contributing to gg — ► J/tpg 
with direct and crossed box diagrams employing the c — c— Jji\> vertex. The crosses indicate that 
the quarks are on-shell. (c) Box diagram with the 4-point vertex c—c—J/iJ) — g. Diagrams with 
reversed quark lines are not shown. 



The partonic differential cross section obtained from the amplitude calculated 
from our model (see [H] for details) is summed over the gluon polarisations, i.e., 

^ = ^— 2 E P^l*. r = L,T 1 ,T 2 , (9) 



where |A^ P9rs | 2 is the squared polarised partonic amplitude for g p {ki)g q {k 2 ) — > 
Q r (P)g s (q) averaged only over colour for polarised cross sections. Here, p, q, r and 
s are the helicities of the respective particles, and s = (ki + k 2 ) 2 , t = {k 2 — q) 2 and 
u = (ki — q) 2 are the Mandelstam variables for the partonic process. The relation 
to the double-differential polarised cross section in transverse momentum Pt and 
rapidity y then is given by 

do r f 1 d ^ 2sP T g(x l )g(x 2 (x 1 )) da r 



dydP T y^mm y/s(y/sxi — E T ey) dt 

In the present calculations, we use the LO gluon distribution g(x) of [21]. 

Figure [4] shows our results with parameter values a = 4 and k = 4.5 GeV for 
a/s = 1.8 TeV in the pseudorapidity range \r)\ < 0.6. The values of a and k were 
chosen to reproduce the cross-section measurement of direct J / ip by CDF (2] up to 
about Pt = 10 GeV. At higher Pt, our curve falls below the data as expected from 
the genuine 1/Pj< scaling of a LO box diagram. We expect higher-order corrections 
incorporating fragmentation- type topologies (~ 1/-Pt) [ES EH] and associated- 
production channels to fill the gap between data and theory at high Pt [26] . 
Figures [H] show our results at a/s = 200 GeV, still with a = 4 and k = 4.5 GeV, 
compared with the PHENIX dats[j] (27] for the central rapidity region \y\ < 0.35 
(a) and the forward one 1.2 < \y\ < 2.2 (b). 

J/ip polarisation measurements at the Tevatron exist only for the prompt yield, 
we have thus computed a from our direct- J/ip cross sections in two extreme cases, 



Note that the PHENIX analysis deals with the total J/tp yield, whereas our computation is for 
the direct yield. For the PHENIX kinematics, the B feeddown can be safely neglected. To what 
concerns the feeddown from Xc it is likely to affect the polarisation observable a (see later), but 
normally much less the Pt dependence. 




FIGURE 4. Comparison between polarised ( 
[with parameters a = 4, k= 4.5 GeV in Eq. (|8|], 
(■v/i= 1.8 TeV, pseudorapidity |r?| < 0.6). 



r r and a£) and unpolarised (cr to t) cross sections 
and CDF experimental data [5] at the Tevatron 




one where the J/^'s from Xc are 100% transverse and another where they are 
100% longitudinal, the first scenario being the more likely one. Figure [6] shows the 
comparison between this computation and the recent results by CDF at y/s = 1.96 
TeV Figures [7] show the polarisation of the direct yield for the central and 
forward rapidity regions at RHIC. At very small Pr, the J/ip is found to be rather 
transversal, a being systematically larger for larger rapidity. 
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FIGURE 6. Prompt J/i/j polarisation: theory vs. CDF data g]. 
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FIGURE 7. (Right) Direct J/ip polarisation as a function of Pt for the RHIC kinematics in 
the rapidity intervals \y\ < 0.35 and 1.2 < \y\ < 2.2. (Left) Direct J/ip polarisation integrated over 
Pt as a function of the rapidity for the RHIC kinematics. 



4. CONCLUSION 



In [E] , we showed that there exist two singularities contributing to the imaginary 
part of the amplitude for gg — > J/ipg. The first can be identified to the CSM 
contribution when the static limit is taken (no relative momentum between quarks). 
The second can be referred to as an s-channel cut and was never considered 
before [13] . 

To deal with such configurations, we have to introduce a four-point function 
c — c—J/ip — g, complementing the information given by the three-point function 
(or Bethe-Salpeter amplitude) c — c— J/tp. Such a four-point function is a priori 
constrained by a low-energy limit (when the emitted gluon is soft) and a scaling 
limit (when the emitted gluon is hard). Given those two physical constraints, we 
constructed a four-point functions exhibiting a dependence only on two parameters, 
which we fixed to reproduce the Tevatron measurements up to mid Pt- We then 
used the latter to compute the cross section for the RHIC kinematics, for which 
we obtained a striking agreement with PHENIX data. This agreement can be 
employed [28] to investigate on the kinematical effects attributable to the final- 
state-gluon emission in studies of shadowing effects on J/ip production in pA, 



dA and AA collisions, in the spirit of the study [221 ED]- Our prediction for 
the polarisation for the prompt J /if) yield at mid Pt at the Tevatron is mostly 
longitudinal. 

In the COM, colour-octet matrix elements account for transitions between a 
coloured heavy-quark pair into a quarkonium by soft unseen gluon emissions 
in the final state. In the present approach, the 4-point function accounts for 
gluon exchanges between the heavy quarks emitting the final-state gluon. As for 
the matrix elements of NRQCD, which are unknown and then fit, we fixed the 
unconstrained parameters of this function in order to reproduce the experimental 
data at ^/i = 1.8 TeV from the CDF collaboration at the Tevatron for Pt < 10 
GeV. 

Contrary to usual results obtained with LO calculations, our approach agrees 
with data down to very low values of the transverse momentum without need of 
resummation of initial-state gluon effects. This feature could be attributed to the 
threshold associated with the cut in the s-channel and should be analysed in more 
details in the future. 

Before drawing further conclusions, several points have to be addressed: Firstly, 
the size of the real part of the amplitude has to be evaluated. When fixing 
the parameter of our four-point function to describe the Tevatron data we have 
implicitly assumed that such a real part was small; this has to verified. Secondly, 
the four-point function we proposed here has to be applied to other regimes of 
production: a similar enhancement by inclusion of the s-channel cut is expected in 
all production processes where the J /if) is associated with a gluon, e.g., photon- 
photon collision at LEP as well as in photo- and lepto-production at HERA. On 
the other hand, other observables insensitive to the COM or the s-channel cut - 
and thus to the ambiguity attached to the description of the four-point function- 
should be studied in the future, especially at the LHC. To conclude, let us mention 
two promising new observables, J / if) production in association with a cc pair [2S] 
and the hadronic activity around the J /if) [31] • 
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